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A STUDY OF IN VISCID FLOW ABOUT AIRFOILS AT HIGH SUPERSONIC SPEEDS 1 


By A. J. Eggehs, Jr., Clarence A. Syvertson, and Samuel Kraus 


SUMMARY 

Steady flow about curved airfoils at high supersonic speeds is 
investigated analytically. With the assumption that air behaves 
as an ideal diatomic gas, it is found that the shock-expansion 
method may be used to predict the flow about curved airfoils up to 
arbitrarily high Mach numbers, provided the flow deflection 
angles are not too close to those corresponding to shock detach- 
ment. This result applies not only to the determination of the 
surface pressure distribution, but also to the determination of the 
whole flow field about an airfoil. Verification of this observa- 
tion is obtained with the aid of the method of characteristics by 
extensive calculations of the pressure gradient and shock-wave 
curvature at the leading edge, and by calculations of the pressure 
distribution on a 10-percent-thick biconvex airfoil at 0° angle of 
attack. 

An approximation to the shock-expansion method for thin 
airfoils at high Mach numbers is also investigated and is found 
to yield pressures in error by less than 10 percent at Mach 
numbers above 3 and flow deflection angles up to 35°. This 
slender-aitfoil method is relatively simple inform and thus may 
prove useful for some engineering purposes. 

Effects of caloric imperfections of air manifest in ■ disturbed 
flow fields at high Mach numbers are investigated, particular 
attention being given to the reduction of the ratio of specific 
heats. So long as this ratio does not decrease appreciably below 
1.3, it is indicated that the shock-expansion method, generalized 
to include the effects of these imperfections, should be substan- 
tially as accurate as for ideal-gas flows. This observation is veri- 
fied with the aid of a generalized shock-expansion method and a 
generalized method of characteristics employed in forms appli- 
cable for local air temperatures up to about 5,000° Rankine. 

The slender-airfoil method is modified to employ an average 
value of the ratio of specific heats for a particular flow field. 
This simplified method has essentially the same accuracy for 
imperfect-gas flows as its counterpart has for ideal-gas flows. 

An approximate flow analysis is made at extremely high Mach 
numbers where it- is indicated that the ratio of specific heats may 
approach close to 1. In this case, it is found that the shock- 
expansion method may be in considerable error; however, .the 
Bvsemann method for the limit of infinite free-stream Mach 
number and specific-heat ratio of 1 appears to apply with 
reasonable accuracy. 

INTRODUCTION 

Small-disturbance potential-flow theories have been em- 
ployed widely, and for the most part successfully, for pre- 


dicting the pressures (and velocities) at the surface of an 
airfoil in steady motion at low supersonic speeds. Thus the 
linear theory of Ackeret (ref. 1) has proven particularly use- 
ful in studying the flow about relatively thin, sharp-nosed 
airfoils at small angles of attack, while the second-order 
theory of Busemann (ref. 2) has found application when 
thicker airfoils at larger angles of attack were under consid- 
eration. At highfree-streamMach numbers the range of appli- 
cability of any potential theory is seriously limited, however, 
due to the production of strong shocks by even the relatively 
small flow deflections caused by thin airfoils. The assump- 
tion of potential flow is invalidated, of course, by the pro- 
nounced entropy rises occurring through these shocks. 

This limitation on potential theories was early recognized 
and led to the adoption (see ref. 3) of what is now commonly 
called the shock-expansion method. The latter method de- 
rives its advantage over potential theories, principally, by 
accounting for the entropy rise through the oblique shock 
emanating from the leading edge of a sharp-nosed airfoil. 
Consequently, so long as the disturbed air behaves essen- 
tially like an ideal gas, and so long as entropy gradients 
normal to the streamlines (due to curvature of the surface) 
do not significantly influence flow at the surface, the shock- 
expansion theory should predict the pressures at the surface 
of an airfoil with good accuracy — it is tacitly assumed, of 
course, that the flow velocity is everywhere supersonic, and 
that the Reynolds number of the flow is sufficiently large to 
minimize viscous effects on surface pressures. 

The departure of the behavior of air from that of an ideal 
gas at the temperatures encountered in flight at high super- 
sonic speeds has been the subject of some investigation in 
the case of flows through oblique shock waves. In reference 
4, the effects of thermal and caloric imperfections on the 
pressure rise across an oblique shock wave were investigated 
at sea-level Mach numbers of 10 and 20, and it was found 
that these effects decreased the rise by less than 5 percent 
for maximum temperatures up to 3,000° R. (corresponding 
to flow deflection angles up to 24°). This decrease was found 
to be due almost entirely to caloric imperfections, or changes 
in vibrational heat capacities of the air passing through the 
shock wave. The changes in temperature and density of the 
air passing through the wave were affected to a considerably 
greater extent. Subsequently, an investigation was carried 
out by Ivey and Cline up to Mach numbers as high as 100 
(ref. 6), using the results for normal shock waves obtained 
by Bethe and Teller considering effects of dissociation (ref. 


i Supersedes NAOA TN 2<H8, "In viscid. Flow About Airfoils at Hi gh Supersonic Speeds" by A. J. Eggers, Jr., and Olarenco A. Syvertson, 1652, and NAOA TN 2729 "An Analysis of 
Suporsonlo Flow in the Region of the Leading Edge of Curved Airfoils, Including Charts for Determining Surface-Pressure Gradient and Shock-Wave Curvature” by Samuel Kraus, 1962. 
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6). As would be expected, the pressures were found to he 
affected to a somewhat greater extent at the higher Mach 
numbers. 

The extent to which flow in the region of the leading edge 
of an airfoil departs from the simple Prandtl-Meyer type 
has also been investigated at high supersonic airspeeds. If 
the surface is curved, for example, to give an exp andin g 
flow downstream of the leading edge, expansion waves from 
the surface will interact with the nose shock wave, thereby 
curving it and yielding a nonisentropic flow field. This 
flow field may be characterized not only by disturbances 
emanating from the surface but also by disturbances reflect- 
ing to some extent from the shock wave back toward the 
surface. The manner in which these phenomena dictate 
shock-wave curvature and surface pressure gradient in 
ideal-gas flows at the leading edge has been treated by 
Crocco (ref. 7) and more recently by Schaefer (ref. 8), Munk 
and Prim (ref. 9), Thomas (ref. 10), and others. In the 
cases considered by Munk and Prim, it was found that sur- 
face pressure gradients were less (in absolute value) than 
those obtained assuming Prandtl-Meyer flow at the higher 
Mach numbers (i. e., Mach numbers greater than about 3) 
although, generally, by no more than about 10 percent. 
Since curved airfoils are likely to be of fundamental interest 
at high flight speeds (see, e. g., ref. 11), the effects of reflected 
disturbances would appear to merit further investigation 
both at the leading edge and as regards their influence on 
the whole flow field. In addition, it would appear desirable 
to consider effects of gaseous imperfections through the field. 

Such an investigation has therefore been undertaken in 
the present report, using the method of characteristics to 
obtain accurately flow fields and as a basis for obtaining the 
more approximate methods of analysis. The method is 
employed in a generalized form which allows caloric imper- 
fections, as well as entropy gradients, in the flow to be con- 
sidered at temperatures up to the order of 5,000° R. — thermal 
imperfections are neglected (see ref. 4). A 10-percent-thick 
biconvex airfoil is treated at Mach numbers from 3.5 to 
infinity, and the results are compared with the predictions 
of the shock-expansion method, including a simplified form 
of the method applicable to slender airfoils at high Mach 
numbers and a generalized form of the method including 
effects of caloric imperfections. In addition, flow' in the 
region of the leading edge of curved airfoils is considered in 
some detail. Values of the surface pressure gradient and 
shock-wave curvature are presented for a wide range of 
Mach numbers and flow deflection angles. 

SYMBOLS 

a local speed of sound 

c chord 

Ci, Ci characteristic coordinates ( Ch positively inclined 
and Ci negatively inclined with respect to the 
local velocity vector) 

Ca section drag coefficient 

Ci section lift coefficient 

c m section moment coefficient (moment taken about 

leading edge) 


c P 

pressure coefficient, ----- 

Cp 

specific heat at constant pressure 

Cv 

specific heat at constant volume 

K 

curvature 

M 

Mach number (ratio of local velocity to local speed 
of sound) 

P 

pressure ratio, ~ 

V 

static pressure 

a 

dynamic pressure 

R 

gas constant 

s, n 

rectangular coordinates (in streamline direction 
and normal to streamline direction, respec- 
tively) 

T 

temperature, °R. 

t 

time 

V 

resultant velocity 

W 

distance measured from leading edge along airfoil 
surface 

x, y 

rectangular coordinates 

ct 

angle of attack, radians unless otherwise specified 

P 

Mach angle, arc sin radians 

y 

£ 

ratio of specific heats, — 

Co 

(Average value of y is y tt .) 

s 

flow deflection angle, radians unless otherwise 
specified 

r 

angle between shock wave and flow direction just 
downstream of shock wave, radians 

e 

molecular vibrational energy constant, °R. (5,600° 
■ R. for air) 

K 

ratio of shock-wave curvature to that givon by tho 
shock-expansion method 

P 

mass density 

a ' 

shock-wave angle, radians 

* 

ratio of surface pressure gradient to that givon by 
the shock-expansion method 

Ol 

ray angle for Prandtl-Meyer flow, radians 

SUBSCRIPTS 

0 *• 

free-stream conditions 

A,B,CP 

conditions at different points in flow field 

i 

ideal-gas quantities 

N 

conditions at the leading edge immediately down- 
stream of the shock wave 

S 

conditions on streamline 

w 

conditions along airfoil surface 

a ' 

conditions along shock wave 


SUPERSCRIPT 

- vector quantities 

DEVELOPMENT OF METHODS OF ANALYSIS 


GENERAL METHODS 

Method of characteristics. — Two-dimensional rotational 
supersonic flows have been treated by numerous authors with 
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tho aid of the method of characteristics, and various adapta- 
tions of the method have been found which are especially 
suited for studying particular types of such flows. In the 
case of steady flows in which atmospheric air does not behave 
as an ideal diatomic gas, a very familiar and simple form of 
the compatibility equations may be employed. To illustrate, 
■consider the Euler equation 

dV , m 

P~dt = — grad# (1) 

the continuity equation 

div (pF)= 0 (2) 


and the equation for the speed of sound (evaluated at con- 
stant entropy) 



(3)' 


Rewriting equations (1) and (2) in the form of partial dif- 
ferential equations and transforming the resulting expressions 
to the characteristic, or C\, Co, coordinate system, there is 
obtained, upon combination with equation (3), the following 
relations for steady flow: 


S2 P . P -4-— V 0 (4) 

pV 2 W, btfJ+Wi ( ) 

and 

cot ( bp bp\ / bS d5 

P V a \dC l + dC s J + {dC l Wj~° (5) 

A simple addition or subtraction of equations (4) and (5) 
then yields the compatibility equations (see, e. g., ref. 12) 


bC~ P V2tim ^c J ( & ) 

and 

wr pV2 tan p (7) 

Now, in reference 4 both caloric and thermal imperfections 
of air were considered, and it was found that the latter 
imperfections 3 have a negligible effect on shock processes in 
atmospheric air. It may easily be shown that this conclusion 
also applies to expansion processes and, for this reason, 
caloric imperfections only are considered in detail in the 
present paper. These imperfections become significant in 
air at temperatures greater than about 800° R. and first 
manifest themselves as changes in the vibrational heat 
capacities with temperature. Thus, the specific heats c p and 
c c and their ratio y for the gas also change. , The equation of 
state remains, however, 


p=pRT 


( 8 ) 


Furthermore, it readily follows from the differential energy 
equation and these expressions that the speed of sound is 
given by the simple relation 

a?=yRT (10) 

Combining equations (8) and (10) and noting that 
sin jS =ajV there is then obtained 


( 11 ) 

sin 1 (3 v ' 

Hence, on combining this equation with equations (6) and 
(7), it is apparent that the familiar compatibility equations 


bp 

— 2 7 p b5 

(12) 

b(7, 

sin 2/3 bCj 

bp 

2 yp b5 

(13) 

b<? 2 

sin 2/3 b(7j 


also hold for the more general type of flow under considera- 
tion. These equations are basic, of course, to two-dimen- 
sional characteristics theory, and, as will be shown later, 
form a convenient starting point for developing sim pler 
theories of two-dimensional supersonic flow. 

In order to apply equations (12) and (13), it is evident 
that the manner in which y and /3-or M are connected to p 
or S must be known.. Relations implicitly connecting these 
variables at temperatures up to the order of 5000° R may be 
readily obtained from the results of reference 4 by simply 
eliminating the terms therein accounting for thermal imper- 
fections. Thus we have as a function of the local static 
temperature and free-stream conditions 


7=7< 


\ i f7«-l\/0V ^ P 

7< )\TJ (e^-l)* 


(14) 


[l+(7r-l) (f) (?^-T) 2 ' J 

and 

2/T 0 \r^ r\ 

7\27L 2 Wi-lV T 0 )^ T 0 \e srr o—l 


(15) 


For isentropic flow along a streamline, the pressure is related 
to the temperature by the expression 


where 


V ACT,) 

p. 4CO 


A CO* 


e srr — 1 
: (• * ,T \ 


n 


or 


(16) 


(17) 


and tho specific heats are still related to the gas constant by 
the expression 

c r —c c =R (9) 

» Thermal Imperfections usually appear In the form of Intennolecular forces and molecular- 
sire effects and may bo accounted for with additional terms In the equation of state. 


If there is a shock wave in the flow, 3 in particular, a nose or 
leading-edge shock, then the folldwing additional relations 
obtained with equations (8), (10), and (15) and the conditions 

» If there are no shock waves, then the subscript <rln-eqoatk>n (10) can, of coarse, be replaced 
with the subscript 0. 
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for continuity of flow and conservation of momentum along 
a stre amlin e through, the shock are also required: 

-rp- (1+y<34 0 2 )+ 

y[(l +y (1 +yoM 0 *) J+4 gj . 

, Ky,TMMy 0 ToM^]-i 
sin (Po/P.y- 1 

and 

tan 6 — tan yo M 0 2 
(P'/Po )— 1 

By use of the local static temperature as a parameter, the 
term 2ypjsm 2)3 in equations (12) and (13) may now he 
evaluated with equations (14) through (17). Equations (18) 
through (20) define the initial conditions downstream of a 
leading edge or other shock wave in the flow field. Thus, 
equations (12) through (20) provide all the information 
necessary to calculate the flow about an airfoil by means of 
the method of characteristics. As described in detail in 
Appendix A, the calculations are of three general types: 
(1) determination of conditions at a point in the flow field 
between the shock and the surface; (2) determination of 
conditions at a point on the surface; and (3) determination 
of condition's at a point just downstream of the shock. Case 
(1) entails the use of both compatibility equations, while 
case (2) entails the use of the compatibility equation for a 
second-family characteristic line in combination with the 
equation of the airfoil surface, and case (3) involves the 
compatibility equation for a first-family line in combination 
with the oblique-shock equations. With the aid of the three 
general types of calculations, the entire flow field about an 
airfoil can be built up numerically using a computing pro- 
cedure working from the leading edge downstream. In cases 
where changes in the vibrational heat capacities with tem- 
perature are neglected, the calculations are, of course, 
simplified since y of the gas can be considered constant, and 
temperature, pressure, and density ratios are simply the 
ideal-gas functions of Mach number. 

Shock-expansion method. — This method of calculating 
supersonic flow of an ideal gas at the surface of an airfoil is 
well known, entailing simply the calculation of flow at the 
nose with the oblique-shock equations and flow downstream 
of the nose with the Prandtl-Meyer equations. Determina- 
tion of airfoil characteristics in this manner requires only a 
small amount of time, of course, compared to that involved 
when the method of characteristics is used, hence, the 
advantage of the former method. The questions arise, 
however, as to exactly what the simplifying assumptions 
underlying the shock-expansion method are, and what form 
the method takes (for calculative purposes) when the gas 
displays varying vibrational heat capacities. 

The matter of simplifying assumptions may perhaps best 
be considered by employing equations (12) and (13), the 
basic compatibility equations. If these expressions are 


(18) 

(19) 

( 20 ) 


resolved into the streamline direction and combined, noting 
that 


and 


bp 1 

bs 2 cos jS 

bS 1 

bs 2 cos jS 



( 21 ) 

( 22 ) 


there is then obtained the relation 

" b8/bC i 

bp bS/bCj 

bs~ 1 . bSfbCi 
[ l+ bt/bO, 

defining the gradient of p along s. If flow along streamlines 
downstream of the nose is of the simple Prandtl-Meyer typo, 
however, we have 


2yp OS 
sin 2)3 bs 


(23) 


Op 2 yp bS 

bs sin 2)3 Os 


(24) 


Hence, it is evident that the requirement for this typo flow is 


bS/bOi 

bSfbGi 


-Cl 


(26) 


Equation (25) is; of course, simply an approximate statement 
of a well-known property of Prandtl-Meyer flows; namely, 
that flow inclination angles are essentially constant along 
first-family Mach lines. It follows from equation (12) that 
if equation (25) holds, then the pressures will also ho essen- 
tially constant along these lines. It' does not follow, how- 
ever, that the. Mach number will be constant or, for that 
matter, that the first-family characteristic lines will be 
straight (as is the case for isentropic expansion flows about 
a corner) . In fact, it may easily he shown that the Mach 
number gradient along C x is proportional to the local entropy 
gradient normal to the streamlines and that the C\ lines are 
curved according to the change in M. Thus we see that 
there is really only one basic assumption underlying the 
shock-expansion method; namely, disturbances incident on 
the nose shock (or, for that- matter, any other shock) are 
consumed almost entirely in changing the direction of the 
shock. In this regard, it is interesting to noto that the 
assumption of Thomas (ref. 13), that pressure is a function 
only of flow deflection angle and entropy, is equivalent to 
this assumption. It follows, of course, that the most 
general solution obtainable with Thomas’ series representa- 
tion of the pressure is that given by the shock-expansion 
method. 

With the assumption that all disturbances incident on the 
shock wave are consumed, it is evident that the shock- 
expansion method provides a relatively simple means for 
calculating the whole flow field about an airfoil, including 
the effects of shock-wave curvature. The details of such 
calculations are presented in Appendix B. In general, of 
course, the validity of this assumption and the accuracy of 
the shock-expansion method can only he checked by com- 
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parison of calculations using this method with those using 
the method of characteristics. 

The shock-expansion method for a calorically imperfect, 
diatomic gas is readily deduced from the equations previ- 
ously obtained. For example, flow conditions at the leading 
edge of an airfoil can be evaluated with the. oblique-shock- 
wave expressions (eqs. (18) through ( 20 )) and the expression 
for conservation of energy (eq. (15)). The variation of flow 
inclination angle with pressure along the surface is then 
obtained by graphically integrating equation (24); namely, 

. „ rPssin20 j 

5s — 5jy — ) — 7)~— &V (26) 

Jpjv p 

where the variables 7 , p, and /S are evaluated using equa- 
tions (14) through (17), employing the static temperature 
as a parameter. When extreme accuracy is not essential, 
this rather tedious calculation can be avoided, and a rela- 
tively simple algebraic solution of the flow downstream of 
the nose can be employed . 4 The details of this solution 
are presented in Appendix C. In the special case of flow 
at high supersonic speeds about slender airfoils, the whole 
calculation becomes particularly simple and warrants special 
attention. 

If it is assumed that the local surface slopes axe small 
compared to 1 and, in addition, that the free-stream Mach 
number is large compared to 1 , it follows that <r and j9 axe 
everywhexe small compaxed to 1 . In this case, equation (24) 
takes on the approximate form 


dp 

d8 



(27) 


Furthermore, if it is assumed that 7 is constant at an average 
value 7 „ for a particular flow field (this assumption appears 
reasonable since, in the temperature range up to 5,000° R, 
the change in 7 is less than about 15 percent), then the Mach 
number and pressure may be related by the simple expression 

M=M n (28) 


Equations (27). and (28) combine to yield the differential 
equation 


— (y.+d 




(£) ,T ‘ <*(£>-<“ ew> 


which readily integrates (between N and S) to the form 

(30) 


now denoting 

^=lM N 8 N =j(M,o N ) (31) 

and 

~~—g(Mo8 N ) ' (32) 

P 0 


* The tabulated results of Noyes (ref. 14) may also prove useful in this case for Mach 
numbers up to 3. 


there is obtained from the oblique-shock equations, simplified 
to conform with this analysis,. 


/(MoSar) — 


M 0 s <r jj- 2 — 1 


and 


where 




g(M 0 8^) 


_ 27aMo 2 try 2 — (7 a — 1) 
7a+l 


MoS„+iJl+( 


7a+l 


Mr 


o5„y 


(33) 


(34) 

(35) 


With equations (30) through (35), the pressures on the sur- 
face of an airfoil may easily be obtained. In terms of pres- 
sure coefficient we have 


*-db[(£XfcH M 

or 

[WCM.W (l-gTp-lj (37) 

The advantage of these slender-airfoil expressions lies, of 
course, in their relative simplicity and, thus, the ease of 
calculation which is inherent to them. It may be noted in 
this regard that the functions / (M 0 5 A .) and g(AI a o M ) can 'be 
calculated once and for all with equations (33), (34), and (35), 
provided the variation of y a with M 0 8 N is known. This 
calculation has been carried out for a constant value of 7 
equal to 1.4 and for average values of 7 , assuming T o =5O0° 
R . 6 The results axe presented in table I. 

It should also be noted that the slender-airfoil expressions 
of the shock-expansion method satisfy the hypersonic simi- 
larity law for airfoils first deduced by Tsien (ref. 15). 8 A 
necessary condition for the validity of these expressions is 
thus satisfied; however, the accuracy of the shock-expansion 
method, whether for slender airfoils or otherwise, remains to 
be investigated. 

METHODS FOR CALCULATING THE FLOW IN THE REGION OF THE 

LEADING EDGE 

As was pointed out previously, pressure disturbances 
emanating from the surface of a curved airfoil interact with, 
and thereby curve, the leading-edge shock wave. The 
geometry of this phenomenon near the leading edge of a 
convex airfoil is illustrated in figure 1 . The pressure dis- 
turbances from the airfoil (expansion waves for a convex 
airfoil) travel along first-family Mach lines Cj. In addition 
to changing the inclination of the shock wave, the interaction 
between these disturbances and the shock produces another 
system of disturbances which travel along second-family 
Mach lines C 3 from the shock wave to the surface. 

Method of characteristics. — An exact solution for the sur- 
face pressure gradient and shock-wave curvature at theleading 


5 For a given valne of To, T.v, to the accuracy of this analysis, is the ideal-gas function 
of Thus, knowing T.V, can bo determined. The average value of y used is 

Y-=T-GUcAy)=('VW-Hv)/2. 

• This fact was employed by Lfnnell (ref. 16) to obtain an expression for pressure coefficient 
equivalent to equation (37) for the caso of constant y and to obtain explicit solutions for the 
lift, drag, and pitching-moment coefficients of several airfoils at hypersonic speeds. 
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Figure 1. — Schematic diagram of supersonic flow past a curved 
sharp-nose airfoil. 


edge may be determined in the following manner. It is clear, 
referring to figure 1, that the change in flow angle between 
points A and C given by the compatibility equations (see eqs. 
(12) and (13)) along the path ABC must equal the change 
determined by the airfoil surface from A to C. ■ Similarly, 
the difference in pressure between points B and D given by 
the compatibility equations along the path BCD must equal 
that determined by the change in shock-wave inclination 
between B and D. In reference 9, these conditions were 
employed at the leading edge to obtain equations, in simple 
parametric form, for determining the surface pressure 
gradients and shock-wave curvatures. These equations can 
be written in the form 


1 dP 
2^ dW~ 


bS/dCP' 

'55/5(7, 

55/5(7, 


2 yP 


sin 2/3 


1 55/5CV 

for the surface pressure gradient, and 


Kv sin (fi ~ a+S ' )+ (M^) 8in ° 5+<r-5) 


(38) 


(39) 


for the shock-wave curvature, where 


55/5(7, 

r 2 yP fdP\ 
_sin 2/3 \d5^ 

~ 

55/5 (7a 

2yP _j_ 
_sm 2/3 

^dSj 

c_ 


T sin (/3 — g-+5) ~l 
Lsin 03+<r— 5)J 


(40) 


It should be pointed out that equation (38) is, of course, 
equivalent to equation (23). 


A procedure for evaluating equations (38) and (39) for a 
calorically imperfect, diatomic gas, as well as for an ideal 
gas, is presented in Appendix D. Since they are exact for 
two-dimensional, steady, inviscid flow, values of the surface 
pressure gradient and shock-wave curvature at the leading 
edge, determined using these equations, may be put to two 
uses. First, the accuracy of approximate methods of cal- 
culating the flow field about a curved airfoil can be evaluated 
at the leading edge by comparing the values of the surface 
pressure gradient and shock-wave curvature predicted by 
these approximate methods to the values obtained using 
equations (38) and (39). Second, the pressure gradient and 
shock-wave curvature can be used to calculate the initial 
points of a characteristic solution for the flow about an airfoil. 

Shook-expansion method. — One approximate solution for 
the surface pressure gradient and shock-wave curvature has 
already been indicated in the previous discussion of the shock- 
expansion method. The requirement for the application of 
the shock-expansion method is .given by equation (25). 
Hence, it is apparent that under this condition the expres- 
sion for the surface pressure gradient (eq. (38)) reduces to 


1 dP 2 yP 

K w dW sin2/3 


(41) 


Similarly, the. expression for the shock-wave curvature (oq. 
(39)) reduces to 


K„_ sin (/3— v+S) 


K w 


(£) 


ao \ . _ 

zrz) sin p 


(42) 


It should be realized that the flow field is determined by 
the basic flow equations in conjunction with the shock wave 
and airfoil surface as boundary conditions. Thus, the addi- 
tional requirement for this shock-expansion method of zero 
pressure gradient along first-family Mach lines means that 
one of the flow relations cannot bo satisfied exactly (i. o., the 
flow field is overdetermined). Equations (41) and (42) 
satisfy the shock relations and the airfoil surface ns boun- 
dary conditions'; however, the compatibility equations are 
only approximately satisfied. (See Appendix B.) 

The error in surface pressure gradient associated with 
neglecting the reflected disturbances might be expected to 
be largest in the region of the leading edge of a curved airfoil 
due to the close proximity of the shock wave and the surface. 
The magnitude of the error in this region may bo deduced, 
of course, from the ratios of values of surface pressure 
gradient and shock-wave curvature given by the character- 
istics method to those given by the shock-expansion method. 
The surface-pressure-gradient ratio and shock-wave-curva- 
ture ratio can be written (using eqs. (38) and (41)) 


4> 


l 


l 


55/5(7, 
55/5 (7 2 
.55/&C7, 
^ 55 / 5(72 


(43) 


and (using eqs. (39) and (42)) 


sin (/3— «r+8)+ (|^^) sin (fi+a—8) 


0 




(44) 
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respectively. A procedure for evaluating equations (43) and 
(44) for the flow of a calorically imperfect gas, as well as for 
an ideal gas, is presented in Appendix D. (The application 
of eq. (43) for an ideal gas has already been given in ref. 9.) 

As was stated previously, the pressure gradient and shock- 
wave curvature can be used to calculate the initial points of 
a characteristic solution for the flow about an airfoil. It is 
apparent that the initial points for a shock-expansion solu- 
tion can be found in a similar manner. With either type of 
solution, flow conditions at an initial point on the surface 
downstream of the leading edge can be calculated with the 
aid of the appropriate value of the surface pressure gradient. 
Similarly, the flow conditions at an initial point on the shock 
wave can be obtained with the aid of the corresponding 
value of the shock-wave curvature. Additional points can 
be obtained between these two by linear interpolation. If 
the two initial points are chosen as the ends of a first-family 
Mach line, there is sufficient information available to deter- 
mine the curvature of this Mach line. Therefore, if detailed 
Icnowledge of the flow in the region of the leading edge is re- 
quired, the surface, shock wave, and first-family Mach line 
can be approximated by circular arcs. (See ref. 17.) 

INVESTIGATION OF THE FLOW ABOUT AIRFOILS AND 
DISCUSSION OF RESULTS 

This study is divided into two sections: first, an investiga- 
tion of the flow in the region of the leading edge of curved 


airfoils; and second, an investigation of the complete flow 
field about an example airfoil. Each of these sections is 
further subdivided into a consideration of the effects of 
Mach number, assuming air behaves as an ideal gas, and 
into a consideration of the combined effects of Mach number 
and gaseous imperfections. In the latter regard, the prin- 
cipal emphasis is placed on the caloric imperfections pre- 
viously discussed. 

FLOW IN THE REGION OF THE LEADING EDGE OF CURVED AIRFOILS 

Ideal-gas flow. — The results of the calculations (using 
eq. (1)4)) of the surface pressure gradient are presented in 
table II and figure 2. 7 The values presented in the table are 
for a range of Mach numbers from 1.5 to y and for leading- 
edge deflection angles from 0° to 45°. Where no value 
appears in the table, the flow behind the shock wave is sub- 
sonic. Corresponding results of the calculations of surface- 
pressure-gradient ratio are presented in table HI and figure 3. 
Erom these results it is seen that except near shock detach- 
ment, the surface-pressure-gradient ratio varies only from 
0.98 to 1.02 for Mach numbers less than 4. Therefore, very 
little error will result from the use of the shock-expansion 
method for the surface pressure gradient at these lower 
Mach numbers. Eor Mach numbers greater than 4 (even 

* Charts were also presented for surface pressure gradient, snrface-presjure-gradlent ratio, 
and shock- wave curvature for Ideal-gas flows In reference 9; however, the results given In the 
present report arc somewhat more ei tensive. Theso results were also presented In cross- 
plotted forms In reference 17. 



Leading-edge deflection angle, S, degrees 

(a) Mo, 1.5 to 6 ' (b) llfo, 6 to 20 

Figure 2. — Variation of surface pressure gradient with leading-edge defleotion angles for various free-etream Mach numbers. 
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Figure 3. — Variation of surface-preesure-gradient ratio with deflection 
angle for various free-stream Mach numbers. 

including a>), 8 the ratio only varies between about 0.9 and 
1.1 (again except near detachment) and, therefore, only small 
error might he expected from the use of the shock-expansion 
method. Near detachment, at the higher Mach numbers, 
the surface-pressure-gradient ratio attains a very large range 
of values and the maximum value increases with Mach num- 
ber. For these conditions, then, the use of the shock-ex- 
pansion method for calculating the surface pressure gradient 
at the leading edge would result in appreciable. error. 

The flow along the surface is isen tropic. Hence, it can be 
shown that >p> the ratio of surface pressure gradient to that 
given by the generalized shock-expansion method, is also the 
velocity-gradient ratio, the Mach number gradient ratio, the 
Mach angle gradient ratio, the density-gradient ratio, and 
the temperature-gradient ratio. Any of these gradients may 
be found, then, by calculating the gradient, using the shock- 
expansion method and applying the appropriate value of 4'- 
This property of the ratio ^ makes it useful in the application 
of the method of characteristics with any of the coordinate 
systems commonly employed in the compatibility equations. 

The results of the calculations of shock-wave curvature 


(eq. (D5)) are presented in table IY and figure 4.° Similarly, 
the results of the calculations of shock-wave-curvaturo ratio 
(eq. (44)) are presented in table Y and figure 6. Except near 
detachment, the curvature ratio varies from 0.92 to 1.08 for 
all Mach numbers including <= . 10 Thus, only small errors 
would result from using the value of the shock-wave curva- 
ture given by the shock-expansion method for all flow con- 
ditions except near shock detachment. 

CaloricaJly-imperfect-gas flows. — With increasing Mach 
number and leading-edge slope, the temperature ratio across 
an oblique shock wave increases as shown in figure 6. As 
the temperature behind the shock wave increases, the be- 
havior of air diverges from that of an. ideal gas as discussed 
previously. Below 800° R., the’ divergence is not significant, 
and the equations for ideal-gas flow can be applied with only 
minute errors resulting. Above 800° E.., the energy of the 
vibrational degrees of freedom of the gas molecules is appre- 
ciable and becomes greater with increasing temperature. 
For these conditions, the specific heats and their ratio vary, 
significantly with temperature. The equations developed 
previously consider these effects and permit the extension of 
the solution for surface pressure gradient and shock-wave 
curvature to the case of calorically imperfect gases. These 



Figure 4. — Variation of shock-wave curvature with loading-cdgo 
deflection angle for various free-stream Mach numbers. 


* For the particular cas e of infinite free-stream Mach number and zero deflection angle, the 
pressure-gradient ratio is double valued- From equation (D6), it is apparent that ^ Is unity 

or zero deflection. Yet, at infinite Mach number. ^ approaches — ^ - v 2 ? — as the 
deflection angle approaches zero. 


• As In the case of K,fKw has the double values of 0 and for Afo« », i-0. 

?4-l 

As in the case of $> and KgfKw, c has the double values of 1 and : 

jliOF" co t 3—0. 


2(2r— 1)[X- V(7'-l)/2rJ 


for 





Leading-edge deflection angle, 8, degrees Shock-wave-curvature ratio, 
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Leadjng-edge deflection angle, 8, degrees 


: 6. — Variation of shock- wave-curvature ratio with leading-edge 
deflection angle for various free-stream Mach numbers. 



Figure 6 . — Variation of leading-edge deflection angle with free-stream 
Mach number for various values of the temperature ratio. 


equations are valid for temperatures up to the order of 5,000° 
R. For a free-stream temperature of 500° R., therefore, the 
shaded area between, lines of constant temperature ratio of 
1.5 and 10, in figure 6, represents the range of conditions for 
which the method developed in this report for the flow of a 
diatomic, calorically imperfect gas would apply. 

The excitation of the vibrational degrees of freedom of the 
gas molecules requires a finite number of collisions, causing 
the well-known heat-capacity lag discussed in references 
5 and 6. The flow distance (i. e., along the streamline) 
required to establish equilibrium is usually small in dense air 
and will be considered infini tesimal in this report. Also, the 
dissociation of air (see ref. 6) will not be considered here. 

Since the free-stream static temperature is an additional 
parameter in calculations of flow of imperfect gases, only a 
limited number of calculations of (1 /K v ) ( dP/dW ), f, Iu/Ku, 
and k were made. The purpose of these calculations is to 
compare the variations of these quantities with the values 
as given by the ideal-gas-flow computations. The calcula- 
tions followed the procedure described in Appendix D. A 
free-stream static temperature of 500° R. was used. The 
results of these calculations are presented in table VI for 
various Mach numbers and leading-edge deflection angles. 

The surface pressure gradients for an ideal gas and for a 
calorically imperfect, diatomic gas are compared in figure 7. 
In all cases calculated, the gradient is smaller for the im- 
perfect gas and diverges gradually, with increasing free- 
stream Mach number and deflection angle, from the value 
of the gradient for an ideal gas. This divergence is consistent 
with the increasing effects of the caloric imperfections due 
to the increasing temperature behind the shock wave. 

The surface-pressure-gradient ratio for the imperfect gas 
is compared in figure 8 with the ratio for an ideal gas. The 
surface-pressure-gradient ratio is smaller for imperfect-gas 
flows than for ideal-gas flows indicating that the effects of 
sliock-wave and expansion-wave interaction are greater. 
This result is attributed, in part, to the fact that the angle 
between the shock wave and the airfoil surface is smaller for 



Leading-edge deflection angle, 8, degrees 


Figure 7. — Comparison of the variation of surface pressure gradient 
with leading-edge deflection angle for various free-stream Mach 
numbers for an ideal gas and for a calorically imperfect, diatomic gas. 
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Leading-edge deflection angle, 8, degrees 


Figure 8. — Comparison of the variation of 6urface-pressure-gradient 
ratio with leading-edge deflection angle for various free-stream Mach 
numbers for an ideal gas and for a calorically imperfect, diatomic gas. 



0 10 20 30 40 50 

Leading-edge deflection angle, 8, degrees 


Figure 9. — Comparison of the variation of the shock-wave curvature 
with the leading-edge deflection angle for various free-stream Mach 
numbers for an ideal gas and for a calorically imperfect, diatomic gas. 



Leading-edge deflection angle, 8, degrees 


Figure 10. — Comparison of the variation of shock-wave-curvature 
ratio with leading-edge deflection angle for various free-stream Mach 
numbers for an ideal gas and for a calorically imperfect, diatomic gas. 


imperfect-gas flows than for ideal-gas flows. The difFeronco 
between the imperfect- and ideal-gas calculations increases 
with increasing Mach number and deflection angle, as might 
be expected. In the range of Mach numbers and flow deflec- 
tion angles investigated, however, the extreme values of ip 
differ only by about 5 percent (see fig. 8). Thus, it is appar- 
ent that while the shock-expansion method •will not bo quite 
as accurate for caloricafly-imperfect-gas flows as for ideal-gas 
flows, the method will not be expected to be invalidated. 

A divergence with Mach number and deflection angle is 
also apparent in figure 9 in which the shock-wave curvatures 
for an ideal gas and a calorically imperfect, diatomic gas are 
compared. This divergence is compatible with the change 
in surface pressure gradient due to the caloric imperfections 
of the gas. The shock-wave-curvature ratio for a calorically 
imperfect, diatomic gas and this ratio for an ideal gas are 
shown in figure 10. Again, it is seen that the effect of caloric 
imperfection is to increase the effects of shock-wave and 
expansion-wave interaction. 

' COMPLETE FLOW FIELDS 

Ideal-gas flows. — The effects of Mach number of primary 
interest here are, of course, those which result from the inter- 
action between the leading-edge (or other) shock wave and 
small disturbances originating on the surface of an airfoil. 
Although the reflected disturbances that are the product of 
this interaction will have the largest effect on the flow near 
the leading edge, their influence on the complete flow field 
about an airfoil also warrants investigation. Further in- 
sight into these effects can be obtained in the region just 
downstream of the shock wave without regard for the shape 
of the airfoil producing the shock. To this end, it is con- 
venient to consider the ratio ~ dp/d^ jjl (see eq. (26)) 

which may be termed “the disturbance strength ratio” since, 
in the region under consideration, it is a measure of the ratio 
of strengths of disturbances reflected from the shock wave 
to disturbances incident on the wave. This ratio may be 
evaluated with equation (40). This calculation has been 
carried out for Mach numbers from 3.5 to <» ( 7 = 1 . 4 ) and 
flow deflection angles approaching those corresponding to 
shock detachment (i. e., M„ <=> 1), and the results are pre- 
sented in figure 11. It is evident that except near M,—\, 



Deflection angle, 8, degrees 


Figure 11. — Variation with deflection angle of the disturbance strength 
ratio behind an oblique shock wave for various free-stream Maoh 
numbers (■y=1.4). 
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the ratio is small (in absolute value) compared to 1 through- 
out the entire range considered — this observation also ap- 
plies, of course, at lower supersonic Mach numbers. Thus 
it is indicated that almost all of on incident disturbance is 
generally absorbed in the shock wave, provided the air be- 
haves like an ideal diatomic gas. This result is substan- 
tially the same, of course, as that which is assumed in deriv- 
ing the shock-expansion method of calculating flows about 
airfoils, and therefore yields additional credence in the method 
for high Mach number os well as low Mach number applica- 
tions. It should also be noted that this result is contrary 
to that obtained by Lighthill (ref. 18), who reports that for 
hypersonic flows(M 0 5^>1) a disturbance is reflected from 

strength. Lightliill’s conclusion is based on an incorrect 
evaluation of his results for the case of very high Mach 
numbers. 

As an over-all check on the shock-expansion method, 
surface pressure distributions calculated thereby are com- 
pared in figure 12 with those obtained with the method of 
characteristics for a 10-percent-thick parabolic-arc biconvex 
airfoil (a=0°) operating at free-stream Mach numbers of 
3.5, 10, and (Additional calculations presented in ref. 
19 were also performed for M 0 — 5, 7.5, and 15.) Predictions 
of the slender-airfoil approximation to the former method 
for high supersonic speeds are also shown. There is no 
apparent difference between the pressure distributions given 
by the method of characteristics and the shock-expansion 
method at a Mach number of 3.5; at a Mach number of 10, 
and more so at infinite Mach number, however, the latter 
method predicts pressures which are slightly low downstream 
of the nose. This result would be deduced from figure 11 
where it is observed that, at the Mach numbers under con- 
sideration, expansion waves incident on the nose shock wave 
are reflected back toward the surface as compression waves 
of relatively small but increasing strength with increasing 
Mach numbor. The effect of these waves does not become 
pronounced even at infinite Mach number (see fig. 12 (c)) 
and it can be seen, upon comparison of these results with 
those presented previously for the pressure-gradient ratio 
(see fig. 3) that the effect of the reflected waves dissipates 
somewhat downstream of the nose. The shock-expansion 
method is thus further substantiated as being a reliable 
simplified method for predicting the flow about airfoils at 
high supersonic speeds, again, so long as the air behayes as 
an ideal diatomic gas. The further simplified slender-airfoil 
method also appears to be a good approximation over the 
entire range of Mach numbers, 11 although, as would be ex- 
pected from the assumptions made in its development, it is 
in somewhat greater error than the shock-expansion method 
at lower Mach'numbers. 

The shape of the shock wave given by the shock-expansion 
method, as presented in Appendix B, is compared in figure 
13 with the shape given by the method of characteristics 
for the biconvex airfoil at M 0 = » . (These shock waves 
correspond to the pressure distributions given in fig. 12 (c) .) 
The shock-expansion method gives a reasonably good 

u The hybrid expression tor pressure coefficient obtained by Ivey and Ollne (re{. 6) gives 
reasonably good results also, although not as accurate as the slender-airfoil method at the 
higher Mach numbers under consideration. 



Chordwise station, x/c 


(a) Mo— 3.5 

(b) Afo=10 

(o) Hfo= » 

Figure 12. — Pressure distribution on 10-percent-tbick biconvex air- 
foil section for various free-stream Mach numbers at a=0°. 

approximation to the shock-wave shape though, as would 
be expected from the results given in figure 5, the curvature 
is somewhat too small. (A procedure for getting a closer 
approximation to the shock-wave shape is also given in 
Appendix B.) It can be seen, however, that this method of 
determining the shock-wave shape is far better than the 
assumption of a straight shock wave that is often associated 
with the shock-expansion method. Evidently, then, the 
shock-expansion method can also be used to calculate the 
flow in regions away from the airfoil surface. (See Appendix 
B.) 
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Figure 13. — Shape of the shock wave for 10-percent-thick biconvex 
airfoil section at <*=0° and M 0 = <*> . 


The relative accuracy at high Mach numbers of the 
slender-airfoil method and linear and second-order potential 
theories may be seen in figure 14. As might be expected, the 
slender-airfoil method is more accurate than linear theory 
at both M 0 —5 and 15 and more accurate than second-order 
theory at M 0 — 15. It is perhaps surprising to note, however, 
that at the lower Mach number of 5, the slender-airfoil 
method is also somewhat superior to the second-order theory. 

The pressure distributions of figures 12 and 14 have been 
employed to calculate the zero-lift drag of the biconvex 
airfoil, and the results of these calculations, along with 
additional predictions of the different methods, are shown 
in figure 15. Predictions of the shock-expansion method are, 

. of course, in best agreement with those of the method of 
characteristics; while the slender-airfoil method, although 
slightly less accurate than the shock-expansion method, is 
apparently superior to both linear and second-order theories 
at Mach numbers above 3.5. 

The preceding findings verify that, so long as the dis- 
turbance strength ratio is small compared to 1, the flow 
along streamlines is essentially of the Prandtl-Meyer type. 
If we choose, on the basis of these findings, a maximum 

absolute value for of 0.06 (note the maximum value of 

O8/OG2 

f or the cases presented in fig. 12 was approximately 

0.06 at M a = <»), the region in which the shock-expansion 
method is applicable can readily be obtained from figure 11. 
The upper boundary line of this region is shown in figure 16, 
and it is evident that it lies only slightly below (about 1°, 
in general) the lin e corresponding to shock detachment given 
approximately by the M, ,=1.0 line. Alm ost the entire 
region of completely supersonic (ideal gas) flow is then 
covered by the method. (See shaded area' of fig. 16.) The 
range of applicability shown in figure 16 is appreciably 
larger than that indicated by Kand (ref. 20) who required 



(a) M 0 =5 

(b) M 0 =15 

Figure 14. — Pressure distribution on 10-percent-thick biconvex airfoil 
section at a=0°. 

that the entire flow field be of the true Prandtl-Meyer type 
(i. e., that all flow properties be constant along first-family 
Mach lines and not just 5 and p). The results presented in 



Figure 15. — Variation of drag coefficient with frec-stream Maoh. 
number for 10-percent-thick biconvex airfoil seption at a=0° (7=1.4) . 
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I 2 3 4 5 6 7 8 9 10 II 12 13 14 15 16 17 18 

Free-slream Mach number, M 0 

Figure 16. — Range of applicability of shock-expansion method 


(7-1.4). 


figures 12 and 13 show, however, that this restriction is not 
necessary. 

The question naturally arises concerning the corresponding 
range of applicability of the slender-airfoil method. This 
question may be answered, in part, by comparing separately 
the predictions of the method for oblique-shock flows and 
expansion flows with those of the exact oblique-shock 
equations and Prandtl-Meyer equations. Such a com- 
parison is shown in figure 17 in terms of the percentage 




(a) <*=10.0° 

(b) <*=19.9° 

(c) <*=30.0° 

Figure 18. — Pressure distribution on 10-percent-thiok biconvex airfoil 
section for various angles of attack at Af 0 =10 (■y=1.4). 

error in the pressure coefficients predicted by the slender- 
airfoil method. As would be expected, this method does 
not exhibit good accuracy over the wide range of appli- 
cability of the shock-expansion method; however, it is 
indicated that it should predict pressure coefficients with 
less than 10-percent error down to Mach numbers as low 
as 3 for airfoils producing flow deflections up to as high as 25°. 

vis a further check on the utility of the slender-airfoil 
method, the pressure coefficients on the 10-percent-thick 
biconvex airfoil have been calculated with this method and 
the shock-expansion method at a Mach number of 10 and 
angles of attack up to about 30°. 12 The results of this 
calculation are shown in figure 18 (see fig. 12 (b) for a=0°) 


(a) Oblique-shock-wave flows. 

(b) Expansion flows. 

Figure 17. — Acouracv of slender-airfoil method in predicting pressure 
coefficients (y=1.4). 


» These candltkms ere within the range of applicability of the shock-expansion method 
as defined In figure 16; hence, the use of the method as a base of comparison seems justified. 
Since the shock-expension method Is far less tedlons to apply than the method of character- 
istics, it will be employed aa such a beso in subsequent calcnlatlans whenever the conditions 
being Investigated have been determined to be within Its range of applicability. 
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(a) Lift coefficient. 

(b) Drag coefficient. 

(c) Moment coefficient. 

Figure 19. — Variation of force and moment coefficients with angle of 
attack for 10-percent-thick biconvex airfoil section at Af o =10 
(7=1-4). 


where it is seen, that the agreement is reasonably good, even 
at the highest angle of attack. This fact is reflected in 
figure 19 showing the force and moment coefficients for the 
airfoil as a function of angle of attack. Little difference 
is observed in the force coefficients as calculated by the 
two methods, while the moment coefficients display more 
pronounced but, nevertheless, small differences at the 
higher angles of attack. 

From these and previous considerations, the ranges of 
applicability of the shock-expansion and slender-airfoil 
methods for supersonic ideal-gas flows are reasonably well 
established. It remains now to determine the manner and 
extent to which gaseous imperfections in the flow at higher 
supersonic speeds may alter these ranges, and the reasons 
therefor. 

Imperfect-gas flows.' — As a first step toward investigating 
the effects of gaseous imperfections on the high Mach 
number flows under consideration, it is convenient to extend 

our consideration of the disturbance strength ratio 

It is recalled that when air exhibits a constant value of 7 
equal to 1.4 (the value for an ideal diatomic gas), the disturb- 
ance str eng th ratio is small at arbitrarily large Mach 'num- 
bers, provided the flow deflection angles are not too close 
to those for shock detachment. One of the most important 
effects of gaseous imperfections is, however, to decrease 7 
of the disturbed air below this value due to the excitation of 


additional degrees of freedom (e. g., vibrational) in the 
molecules at the high temperatures encountered at high 
Mach numbers. Indeed, at arbitrarily high Mach numbers 
it might be expected that 7 of the disturbed ah would 
approach 1 , since the numbejr of degrees of freedom may 
effectively become very large (see, e. g., refs. 3 and 6 ). In 
this case, however, , the extent of the disturbance flow field 
is decreased -to a layer at the surface of the body which is 
negligibly thin compared to that for the case of ideal-gas 
flow. Thus, it is apparent that significant changes in the 
flow about airfoils at high Mach numbers may result from 
decreases in 7 of the disturbed air; hence, the effects of such 
decreases on the disturbance strength ratio would appear to 
warrant attention. 

A detailed analysis of these effects is impractical at the 
present time, due to the limited range over which the vari- 
ation of 7 with temperature is accurately known. Even in 
the range where this variation is so known, the additional 
complication required to consider the effects of variable 7 
and the addition of another independent parameter (free- 
stream temperature) make extensive calculations of the 
disturbance strength ratio impractical. However, some 
knowledge of these effects can be gained by performing the 
calculations for one free-stream Mach number and tem- 
perature. Such calculations have been carried out at a 
Mach number of 10 for a free-stream temperature of 600° 
R. and the results are presented in figure 20 . The curve for a 
calorically imperfect gas cannot be extended to Bhock 
detachment because the temperature behind the shock wave 
exceeds that for which the caloricafly-imperfect-gas equations 
are valid. It can be seen that the effect of the caloric im- 
perfection of air is to increase the value of the disturbance 
strength ratio and that the effect increases with increasing 
temperature or decreasing 7 . However, it appears that if 7 
does not decrease appreciably below 1.3, as in this case, the 
disturbance strength ratio is still small compared. to unity. 
It might be expected, therefore (as previously found for flow 
in the region of the leading edge), that the shock-expansion 
method would continue to predict the flow about complete 
airfoils with reasonable accuracy. This point has been 
checked with the methods developed previously for analyzing 
the flow of a calorically imperfect, diatomic gas at local air 
temperatures up to about 5,000° R. (note 7 has a value only 
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Figure 20. — Effect of the caloric imperfections of air on the disturbance 
strength ratio at Af 0 =10 and 3o=500° R.‘ 
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Figure 21. — Pressure distribution on lower surface of 10-percent-thick 
biconvex airfoil section at Mo=10, 2o=500° R, a=19.9°. 


slightly less than 1 .3 at this temperature) . In particular, the 
pressure distribution on the lower surface of the biconvex 
airfoil at M 0 =10, a=19.9°, and T 0 = 500° R. (2V«4000° R.) 
lias been calculated with both the method of characteristics 
and the shock-expansion method. 13 The results of these 
calculations are presented in figure 2 1 , and it would appear that 
the conclusions drawn from figure 20 are substantiated. (Due 
to the lower temperatures, pressures in the expansion flow 
about the upper surface are not influenced by caloric imper- 
fections and, hence, are the same as shown in figure 18 (b).) 
For these same conditions, the shapes of the shock waves 
given by the shock-expansion method and by the method of 
characteristics (both for a calorically imperfect gas) are com- 
pared in figure 22. Just as in the case of ideal-gas flows, the 
shock-expansion method gives a good approximation to the 
shock-wave shape, far better than the assumption of a 
straight shock wave. Thus, it is seen that the conclusion 
drawn from figure 20 should also apply for the use of the shock- 
expansion method to calculate the flow field away from the 
airfoil surface. 

Shown also in figure 21 is the pressure distribution ob- 
tained by the shock-oxpa'nsion method for an ideal gas 
(y 1 =1.4). It is apparent, on comparing this pressure distri- 
bution with the other distributions, that although the effect 
of caloric imperfections on the disturbance strength ratio 
is small, the pressures are appreciably reduced by the increase 
in specific heats. The extent of this reduction is more com- 
pletely illustrated in figure 23 where the lower-surface pres- 
sure distributions on the biconvex airfoil are presented for 
Mo— 10 and T o =500 0 R., at a=0°, 10°, 19.9°, and 30°. As 
one might expect, the reduction in pressures increases with 
angle of attack (due to the corresponding increase in static 
temperature of the disturbed air). The pressure coefficients 
calculated with consideration of the imperfections in the gas 
are less on the lower surface (up to 6 percent at the leading 
edge and 15 percent at the trailing edge) than those calculated 
assuming the gas behaves ideally. The upper-surface pres- 


u For added ease of calculation the expansion mothod of Appendix O was employed. Thl3 
method Is also employed In all subsequent calculations of this type since It has been found to 
yield results dlflerlng by less than 1 percent from those obtained by the more tedious graphical 
Integration method. 


sures are again unaffected by the caloric imperfections of air 
in all the cases presented (except at a=0°) since this surface 
experiences lower pressures and, hence, lower temperatures. 
They are therefore the same as shown in figure 18. Shown 
also in figure 23 are the pressure distributions calculated 
with the slender-airfoil method for y=y a . The accuracy of 
this simplified method is substantially the same as was pre- 
viously observed for the corresponding method in the case 
of ideal-gas flows, although the local error may be greater 
than the reduction in pressure coefficients due to the caloric 
imperfections of air. 

The force and moment coefficients, corresponding to the' 
lower-surface pressure distributions shown in figure 23 and 
the upper-surface distributions of figure 18, are presented in 
figure 24. The reduction in the lower-surface pressures 
leads, of course, to a general reduction in all three coefficients 
(up to about 10 percent for a=30°). The slender-airfoil 
method again predicts these coefficients with surprising 
accuracy. 

In order to further assess the accuracy of the slender-airfoil 
method, some additional calculations were carried out for the 
biconvex airfoil at a=0° and M 0 =20 and 30. The pressure 
distributions for these cases were calculated by the shock- 
expansion method, slender-airfoil method (y=y<d, and 
slender-airfoil method (7=71). These results are presented 
in figure 2&, and it is observed that the use of y a rather than 
7< improves the accuracy of the slender-airfoil method. 
The extent of this improvement in the case of drag coefficient 
is shown in figure 26; it would appear that predictions of the 
slender-airfoil method (7=70) and shock-expansion method 
are in as good agreement as for ideal-gas flows (see fig. 15). 
On the basis of these and previous results, it may be con- 
cluded, then, that not only does the shock-expansion method 
retain its range of applicability when air exhibits caloric 
imperfections, provided 7 of the disturbed air is not appre- 
ciably less than 1.3, but also the slender-airfoil method 
(7=7o) retains its range of applicability. 

It would be surprising indeed, however, if this conclusion . 
continued to apply as 7 of the disturbed fluid approached 1 



Figure 22. — Shape of the shock wave for 10-percent-thick biconvex 
airfoil section at <*=19.9°, Af 0 =10, To=500° R. 
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' Figure 23. — Pressure distribution on lower surface of 10-percent-thick biconvex airfoil section for various angles of attack at iWo= 10 and 70=500° It. 


since, as indicated previously, ™ creases decreas- 

ing y. Although tlie manner in which, y varies with tempera- 
ture is not known in this range, some knowledge of these 
effects can be gained by repeating the ideal-gas calculations 
for constant values of y between 1.4 and 1.0. w Such calcula- 
tions have been carried out for infinite Mach number since, 

in this case, has its maximum value for a given y, 

and the results are presented in figure 27. It is seen that 
except near detachment, the disturbance strength ratio in- 
creases with decreasing y. This increase is slow at first; 

for example, the value of q I s 8 tiH less than 0.1 at 

7=1.3. This result is in agreement with the previous con- 

it since the enthalpy la negligibly small compared to the mass kinetic energy of the undis- 
turbed fluid at the high Mach numbers of interest and, hence, y of this fluid does not Influence 
the flow, this approach corresponds to employing an average value of y for the disturbed 
fluid. 


elusion regarding cases where y is greater than 1.3. How- 
ever, continues to increase as y decreases, and, in 

fact, approaches 1 as y approaches 1. The effect on pressure 
distributions of this increase in the strength of the reflected 
disturbances may also be investigated by using the ideal-gas 
relationships in combination with appropriate values of y. 

The limiting case of infinite free-stream Mach number and 
7=1.0 (for the disturbed fluid, see footnote 14) has already 
been investigated by Busemann (ref. 21) and more recently 
by Ivey, Klunker, and Bowen (ref. 22). In this case, as 
pointed out previously, the shock wave emanating from the 
leading edge remains attached to the surface downstream of 
the leading edge, (this is easily verified with the oblique- 
shock-wave equations), and the disturbance flow field is 
confined to an infinitesimally thin layer adjacent to the 
surface. In addition, the velocity along a streamline down- 
stream of the shock is constant, as may easily be shown with 



A STUDY OF rcmSCTD FLOW ABOUT AIRFOILS AT HIGH SUPERSONIC SPEEDS 


' 355 


.6 


S .4 


.2 


.3 

C 

01 

’o 

£ .2 
CJ 

o 


•--.I 

*» 


-.3, 






















(a) -- 






















/ 









/ 

✓ 









/ 





(b) 


- 














’ " 








a 

(c) 

■ 

■ 

■ 

■ 

a 

■ 

■ 

a 

a 

Shock-expansion method 

Slendei — airfoil method 

■ 

■ 

a 

a 

a 


16 24 

Angle of attack, a, degrees 


32 


40; 


(a) Lift coefficient. 

(b) Drag coefficient. 

(c) Moment coefficient. 

Figure 24. — Variation of force and moment coefficients with angle of 
attack for 10-percent-thick biconvex airfoil section at M 0 =lO, 
ff’o^SOO 0 R. 


the compatibility equations. Surface pressures therefore 
become a simple function of airfoil geometry 

O v =2 sin a 5s+2 cos 5 S J sin S s dx (45) 


7 approaches 1 . For 7 = 1.0 and M a = « the shock-expansion 
method, in turn, predicts a discontinuous pressure distribu- 
tion with a pressure coefficient equal to that of the Busem ann 
theory at the leading edge hut a pressure coefficient of zero 
at all points downstream of the leading edge. Hence, it 
may be concluded that when the free-stream Mach number 
approaches infinity and 7 approaches 1 , the Busemann 
method rather than the shock-expansion method for calcu- 
lating the flow about airfoils should be employed. 



varying, to a first approximation, directly with the square 
of tho component of free-stream velocity normal to the sur- 
face (i. e., the flow is approximately of the Newtonian 
corpuscular type). With this theory, then, and the method 
of characteristics, we can ( get an idea of both the extent to 
which changes of 7 from 1.4 toward 1 will alter surface 
pressures, and the accuracy with which the shock-expansion 
theory predicts the alterations. To this end, figure 28 is 
presented showing the pressure distributions about the 
biconvex airfoil at M 0 — a> as calculated by the several 
methods for different values of 7 . It is observed that, 
whereas the shock-expansion method agrees very closely 
with the method of characteristics for 7 = 1 . 4 , there is a 
large difference at 7=1.05. This, of course, is precisely 
what one would expect from the previous discussion of the 
disturbance strength ratio. On the other hand, if the two 
characteristics solutions and the Busemann method are 
considered in order of decreasing 7 , it is indicated that the 
characteristics solutions approach the Busemann theory as 


(a) M 0 =20 

(b) M 0 =30 

Figure 25. — Pressure distribution on 10-percent-tbick biconvex airfoil 
section at cc=0° and T o =500° R. 



Free-stream Mach number, M a 

Figure 26. — Variation of drag coefficient with Mach number for 
10-percent-thick biconvex airfoil section for a=0° and 7o==500 0 R. 
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Deflection angle, 8, degrees 

Figtjbb 27. — Variation of disturbance strength ratio with deflection 

angle at infinite free-stream Mach number for various values of y. 

CONCLUSIONS 

Inviscid flow about curved airfoils at higb supersonic 
speeds was investigated analytically, first assuming air be- 
haves as an ideal gas, and then assuming it behaves as a 
thermally perfect, calorically imperfect gas. This study has 
led to the following conclusions: 

1. So long as air behaves as an ideal gas, the shock- 
expansion method may be used with good accuracy to pre- 
dict the flow about a curved airfoil up to arbitrarily high 
Mach numbers, provided the flow deflection angles are 
about 1° or more below those corresponding to shock detach- 
ment. This conclusion applies not only to the determina- 
tion of surface pressure distributions, but also to the deter- 
mination of the whole flow field about the airfoil. 

2. An approximation to the shock-expansion method, 
applicable to ideal-gas flows about slender airfoils at high 



(a) 7 = 1.4 

(b) 7=1.05 

Figure 28. — Pressure distribution on 10-percent-thlok biconvex airfoil 
section at M a — co and a=0°. 

Mach numbers, predicts pressure coefficients in error by loss 
than 10 percent for Mach numbers above 3 and flow deflec- 
tion angles up to 25°. ' 

3. So long as caloric imperfections of air do not decrease 
the ratio of specific heats appreciably below 1.3 (correspond- 
ing to air temperatures up to the order of 5,000° R.), the 
shock-expansion method, generalized to include the effects 
of these imperfections, is substantially as accurate as for 
ideal-gas flows. The principal effect of caloric imperfections 
is to reduce pressure coefficients by as much as 15 percent. 
The slender-airfoil method can also be made as accurate as 
for ideal-gas flows by employing an average value of the 
ratio of specific heats. 

4. If the ratio of specific heats approaches 1, as it may at 
extremely high Mach numbers, the shock-expansion method 
can be in considerable error. In this case, the Busemann 
method for flow in the limit of infinite Mach number and 
specific-heat ratio of 1 applies with reasonable accuracy. 


Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif., January 9, 1952. 





APPENDIX A 


METHOD OF CHARACTERISTICS FOR TWO-DIMENSIONAL 
FLOW OF A CALORICALLY IMPERFECT GAS 

In the application of the method of characteristics for a 
calorically imperfect, diatomic gas to the particular problem 
of analyzing the flow about curved two-dimensional airfoils, 
many of the calculations are identical to those encountered 
in the solution of any problem where characteristics theory 
is employed. Since the details of these calculations are well 
known and well reported (see, e. g., ref. 12), they will not be 
repeated here. 

A lattice-point system with an initial-value, numerical 
computing procedure will be used. The form of the com- 
patibility equations to be employed was developed pre- 
viously; 11 however, it is convenient for purposes of calcula- 
tion to substitute the pressure ratio, p/q 0 , into these equations 
and to rewrite them as difference equations. Equations 
(12) and (13) are thus reduced to the following forms: 

and 
where 


It is also convenient to employ several reference curves. 
Thes.e curves can be divided into two groups. The general 
reference curves consist of 7 and A (T) as a function of 
temperature, T. Equations (14) and (17) are used to deter- 
mine these curves. A second set of shock-wave reference 
curves consisting of p/q 0> <r, and 5 as a function of tempera- 
ture, T, are determined by use of equations (18) through 
(20) — the values of T 0 and M 0 are presumed known. 

In the computations three types of points are encountered. 
These are (1) a point in the flow field between the shock 
wave and the airfoil surface, (2) a point on the airfoil sur- 
face, and (3) a point just downstream of the shock wave. 
Each one of these types of points requires a slightly different 
computing procedure and they will be* considered in order. 

POINT IN THE FLOW FIELD BETWEEN THE SHOCK WAVE AND THE 
AIRFOIL SURFACE 

Figure 29 (a) shows a schematic diagram of the system of 
points to be considered in these calculations. Point C is 
the unknown point at the intersection of the first-family 
characteristic line passing through point A and the second- 
family characteristic line passing through point B. Six 
quantities are known at both points A and B, and the prob- 
lem is to calculate these same quantities at point C. These 

>• This form of the compatibility equations (In p and s coordinates) was also used In ob- 
taining some of the characteristics solutions lor Ideal-gas flows. The majority of these solu- 
tions were carried out, however, with the compatibility equations In 0, l, and entropy coordi- 
nates, since It was found that greater accuracy was usually obtained for a given net sire. 
In general, the net sire employed yielded pressures at from 30 to 3« surface polntb on an airfoil 
with a maximum error In tho corresponding pressure coeflldents equal to less than 1 percent of 
the pressure coeflldent at tho leading edge. 


O/ffo )a— (p/2o) A= — \a (So— 8a) 

(Al) 

(pl 2o) c — (pl ?o)n = Xa (8c — S B ) 

(A2) 

_ 27 (p/go) 
sin 2/5 

(A3) 


quantities are x, y, 8, p/q 0) T, and 'T a . The first five quantities 
are of obvious significance. The sixth, is defined as the 
static temperature, just downstream of the shock wave, on 
the streamline passing through the point C. 

The physical coordinates of the point C(x c , y c ) may be 
determined by standard procedures such as those given in 
reference 12. In order to determine the quantity 8 C , it is 
necessary to solve equations (Al) and (A2) simultaneously, 
thus, 

Sc - \I+T a (A4) 

Equation (Al) or (A2) is then used to obtain (p/q 0 )a ■ 

There remains only the problem of determining T c and 
T, c at point C. The temperature is obviously constant 





(a) Point in field. 

(b) Point on surface. 

(c) Point on shook wave. 

Figure 29. — Diagram of point system in the method of characteristics 
* for the two-dimensional flow of a calorically imperfect gas. 
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along the streamline through C. This quantity may there- 
fore be calculated in 'the same manner as the entropy is 
calculated in similar flow fields for ideal-gas processes (see, 
o. g., ref. 12). Furthermore, since the flow along stream- 
lines downstream of the shock wave is isentropic, equation 
(16) may be applied in the form 

( P/Qq)c AgV (^5) 

(2>/2o). c A (To) 

The pressure, (p/q_o)« G , is defined in a manner analogous to 
T 9 , and may thus be determined using the shock-wave 
reference curves and the known value of 2V C . Similarly, 
A(T<r c ) may be determined from the general reference curves. 
The only unknown in equation (A5) then is A {T c ) which 
may now be calculated. Once A (T c ) is determined, T c may 
be determined by again using the general reference curves. 
All six quantities, x c , y c , S c , (pllo)c, T c , and 2V C have now 
been determined. 

POINT ON .THE AIRFOIL SURFACE 

Figure 29(b) shows a schematic diagram of the points to 
be considered in these calculations. The physical coordi- 
nates of point C(xc, He) are first calculated by solving simul- 
taneously the equation of the second-family Mach line 
passing through point B and the equation of the airfoil 
surface. When x c and y c have been determined, S c is readily 
obtained from the equation of the airfoil surface. Equation 
(A2) is then applied to determine (p/g 0 )c- 

Since the airfoil surface is a streamline, T, c is constant 
along the surface and may be evaluated at the leading edge. 
The temperature, T c , may then be determined using equation 
(A5) and the previously described procedure. All six quanti- 
ties, Xc, y c , (jplZo)c, S 0 , T c , and T, , c , are thus determined. 

In the special case of the first point on the airfoil surface 
downstream of the leading edge, the pressure ratio is calcu- 
lated using the pressure gradient evaluated at the leading 
edge. (See section Methods for Calculating the Flow in 


the Region of the Leading Edge, and Appendix D.) Addi- 
tional initial points in this region may be calculated by the 
procedure previously described. 

POINT ON THE SHOCK WAVE 


Figure 29 (c) shows a schematic diagram of the points to 
be considered in these calculations. The physical coordi- 
nates of point C(x c , y c ) are first calculated by solving simul- 
taneously the equation of the first-family Mach line passing 
through point A and the equation of the shock wave linear- 
ized at point D, the last known point on the wave. The 
variation of p/g 0 with 6 along the shock wave may be approxi- 
mated by the relation 


(p/ffo)c-(p/go)i > =[^^2 {5o ~ 8d) (A6) 

In this equation * s ^ ie ra ^® °f change of p/qa 


with 5 along the downstream side of the shock wave evalu- 
ated at point D. Because of the complicated nature of the 
shock-wave equations, it is generally efisier to evaluate 

graphically or numerically from the shock-wave 

reference curves; however, this derivative may also be eval- 
uated from the equations given in Appendix D. Equations 
(Al) and (A6) are solved simultaneously for S a , thus, 



i i 

•S!« 

i i 

S D + (p/qo)A— (p/qo)D 

Ai+ 

~d(p/qo)~ 
L dS 

9 


(AT) 


when S c has been calculated, T c and, in turn, (p/qo)o may be 
determined from the shock-wave reference curves. Since 
point C in this case is just downstream of the shock wave, 
To and T, c are identical. The.six quantities, x c , yc, (p/qo)c, 
5 a , T Cl and T„ c have now been determined. 



APPENDIX B 


SHOCK-EXPANSION METHOD FOR CALCULATING THE FLOW 
FIELD ABOUT AN AIRFOIL 

An initial-value procedure which, is s imil ar to, although 
markedly simpler than, that associated with the method of 
characteristics may be employed to carry out this calcula- 
tion. 16 To illustrate, consider figure 30. With the oblique- 
shock-wave and expansion equations, all fluid properties 
at points A, B, D, and so forth, on the airfoil surface'may be 



FiauRE 30. — Schematic diagram of shock-expansion method for cal- 
culating the flow field about an airfoiL 


calculated in the usual manner. If point B is chosen close 
to point A, the first-family ((70 Mach line connecting B to 
point C on the shock wave may be considered straight and 
inclined at an angle to the free stream equal to jSs+Sb. 
Similarly, the segment of the shock wave AC may be con- 
sidered straight and inclined to the free-stream direction at 
an angle of <r x . Thus, the physical coordinates of point C 
may be easily calculated. Since 5 is assumed constant along 
first-family Mach lines, 8 a is equal to S B . All fluid properties 
at point C may be calculated from this known value of 8 0 
with the oblique-shock-wave equations. In a similar 
manner, the segment D E of the next first-family Mach line 
is considered straight and inclined at f3 D +8 D , and the 
streamline joining points C and E is considered straight and 
inclined at So- The physical coordinates of point E are 
therefore easily obtained. Since the flow along streamlines 
downstream of the shock wave is isentropic and since the 
pressure is also assumed constant along the first-family 
Mach lines (i. e., P e =Pd), the fluid properties at point E 
are readily obtained from the known properties at point C, 
using the isentropic flow relationships. The construction 
of the remainder of the flow field follows in a similar manner. 

As was discussed previously, the assumption that 8 is 
constant along first-family Mach lines is an additional con- 


dition which, in general, overdetermines the flow field. It is 
possible, therefore, to calculate two values of the shock-wave 
angle at each point on the shock, one assuming S is constant 
and one assuming the pressure is constant. These two values 
will differ slightly and it can be shown that they will bracket 
the correct value that would be given by the method of 
characteristics. It is also relatively easy to show that if 
the change in 8 (or p) along C\ given by the corresponding 
characteristic solution is small, 17 the true value of the shock- 
wave angle lies just midway between the two values given by 
the different assumptions. It is apparent then that a 
closer approximation to the shock-wave shape can be easily 
obtained by simply averaging the values of shock-wave 
angle determined by ass umin g 5 is constant and by ass uming p 
is constant. The increase in accuracy is illustrated in figure 3 1 
for the biconvex airfoil at M 0 = <*> . The averaging procedure 
gives a shock-wave shape that is closer to that given by the 
method of characteristics by 60 to 80 percent. The increase 
in computation time (at least for ideal-gas flows) is negligi- 
ble. In this regard, it is interesting to note that the shock- 
expansion solutions require less than 20 percent the computa- 
tion time of the characteristic solutions. 

The shock-expansion method is applicable to the determi- 
nation of the flow not only in the region adjacent to the 
airfoil (whether concave or convex), but also in the region 
downstream of the airfoil; hence, it may, for example, prove 
useful in downwash studies and the like. 



» It Is clear, of course, that an average-value procedure could also be employed. In fact, 
for tbo shock-expansion method, an average-value procedure requires very little additional 
computation since the fluid properties for a system of points can he obtained Independently 
of their physical coordinates. (The difference In the average- or initial-value procedures 
appears only In computing the physical coordinates.) Thus, the slopes at both ends of a 
lino segment are known before the line Is added to the construction. No Iteration, as re- 
quired In characteristic solutions. Is necessary In this case. 


Figure 31. — Shape of the shock wave for 10-percent-thiek biconvex 
airfoil section at cr=0° and Af 0 = <x> . 


u Thenecessaryassumptlon here Is that and 

equivalent to the condition that Is much less than one. 


are nearly equal. 
(Bee eq. (40).) 


This Is 
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APPENDIX C 


APPROXIMATE SOLUTION FOR PRANDTL-MEYER FLOW OF 
A CALORICALLY IMPERFECT GAS 

The following solution is obtained with an analysis similar 
to that used in Meyer’s original paper (ref. 23) . A schematic 
diagram of the subject flow field is shown in figure 32. It is 



Figure 32. — Schematic diagram of Prandtl-Meyer flow around a comer. 


evident that the change in flow-inclination angle for Prandtl- 
Meyer flow can be written as follows: 


8 N — 8— (fitr — — “) (Cl) 


Since the flow is isentropic, a given value of the local pressure 
will determine the Mach angle, /3. The problem, then, is to 
evaluate the angle, w. To this end, the velocity components 
tangential and normal to the first-family Mach lines may be 
expressed in the usual manner in terms of a potential p, thus, 



(C2) 


Equation (15), for conservation of energy, may be written 
in terms of the local temperature as follows: 

F*+ 2 RT+2RT (08) 

The constant, A, is evaluated at the conditions existing up- 
stream of the expansion region; namely, 

^■=w(^ I )sr»+2B2’„( ? ^£ 7 ) (09) 

Equations (C7) and (08) are then combined to yield 

&+$ a ’=-2RT (-F^+J^+A 3 (CIO) 


It was shown previously, however, that 

a 3 =yRT 


(Oil) 


Equations (C6), (CIO), and (Oil) may therefore be combined 
to obtain the following relationship: 


•*+♦-* [l (012) 


or 


( y<+ 1 


(Tt— i Yr 




TC 


Y< 


(013) 


From the imperfect-gas relationship for y We have 


1+(t<— 1) 

7i 

(- 

\Ty 

{ e* ,T 

1 (e i/T - i)» 

7 , , Y<— 1 / 

kTJ 

5 gS/T 

!+ yt \ 

(e ,rr — l) 2 


(014) 


By substitution of this relation into equation (013) there is 
then obtained 


1 dp 

(L — 

r 0 co 


(03) 


It is clear, however, that these components are functions of 
« only; hence, it is convenient to define a new velocity poten- 
tial which is a function of « alone. Such a potential is 


*(«)=f (04) 

The velocity components may then be written in terms of 
this new potential. 

, «=$ (05) 

a=$ a (06) 

The resultant velocity is given by the expression 

(07) 


where 

f ( t ) 


Now 




e/T 

gS/T —1 


i+(t<— i) 
1 + 


(g/T)e elT r, , e/T 

e *rr—\ 

T< -1 ( d/T) 3 e e/T 
7< (e >/T — l) a 



yT_ a? 

7,0 7 ,f?0 


(016) 


(016) 


(017) 


For every value of T/0 there is thus a particiflar valuo of 
a?/y iRB. The function F(d/T) is therefore uniquely deter- 
mined for any value of a 1 since y iRd is, of course, a constant. 
With this point in mind, let 

F(B/D=G(a 3 ly i R8 ) (018) 
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Figure 33. — Accuracy of approximation used in Appendix C to obtain 
solution for Prandtl-Meyer flow of a calorically imperfect gas. 


Figure 33 shows G(a 2 /viR0) plotted as a function of a a /yiR6. 
This curve is approximated with the following simple 
relation: 

<?Wr,ffi»=0.38 -^+0.71-^ (019) 

for0 - 18 <S< 10 

and 


G(a 2 /y t R6 )= 0 (C20) 


for 0< 


y {RO 


<0.18 


/ 

Equation (C19) is also plotted in figure 33 to show the accu- 
racy of this approximation. Consider first the case when G 
is given by equation (Cl9) which is written in the form 


G{a?ly t l te)=Za*+n+ v ^ (C21) 

U/ 

where, obviously, 

£=0.38 /y t Rd \ 

p=0.71 > (C22) 

V — — 0.14 (y t R6)) 

Equation (C21) is substituted into equation (Cl5) and, with 
equations (Ol8) and (C6), the following expression results: 

— -A 2 =0 (C23) 

y t \yt — 1 yiJ t< 

In order to simplify this equation the following substitutions 
are made : 


A.+i + W,s U , 2 x 

(C24) 

1 

t- 

V 


VI 1 
1 

« 

"'I 

col 


sin 2 v ^ V ^ / 

(C25) 



r-V 

(C26) 


and 


4V 


2 _D i cos* r(r„)* 

A{ 

7i 


(C27) 


Equation (C23) then reduces to 


4 cos v 


4 sin 2 r 4 sin 2 v 


cos 4 t(t cos 2 T ( Tu )*+^i-A-^-A=0 (C28) 


This equation is solved for r«, thus, 

— - — (cos t— cos v) 1 ' 2 
V D cos r v ' 


or 


du— 


D cos 7 dr 
2 (cos r— cos v) 1/2 


(C29) 

(C30) 


This expression is readily integrated to obtain the following 
equation relating w to the local velocity: 


a>-co„=jD {2 [E(Jc,z)~E(k,z K )\-[F(k,z)-F(Jc,z„j ] } (C31) 

where 

E elliptic integral of the second kind 
F elliptic integral of the first kind 

t y 

k sin ^ (modulus) 

s . 9in " 1 (sTT© (ampUtude) . 

The procedure for calculating corresponding values of the 
pressure, p, and the deflection angle, S, is straightforward 
with the aid of the preceding equations and may be sum- 
marized as follows: 

1. Calculate A 2 (eq..(C9)) 

2. Calculate D 2 (eq. (C24)) 

3. Calculate v (eq. (C25)) 

4. Assume a value of T, less than T N 

5. Calculate p (eqs. (16) and (17)) 

6. Calculate V 2 (eq. (C8)) 

7. Calculate y (eq. (C14)) 

8. Calculate a 2 (or <K) (eq. (Cll)) 

9. Calculate M and, in turn, fi from V and a 

10. Calculate u 2 (or 4 s (eq. (C7)) 

11. Calculate r (eq. (C26)) 

12. Calculate u (eq. (C31)) 

13. Calculate S (eq. (Cl)) 

This procedure is followed so long as the quantity a 2 /ytRd 
is greater than 0.18. (This is equivalent to the temperature 
being greater than approximately 1,000° R.) For values 
of a 2 ly' t RQ less than or equal to 0.18 (or temperatures less 
than about 1,000° R.), (?is set equal to zero (see eq. (C20)) 
In this case, equation (C15) reduces to the same form as for 
an ideal gas and, therefore, the well-known ideal-gas rela- 
tionships can be used. 




APPENDIX D 


EVALUATION OF (l/K w )(dP/dW), KJKw, AND * 

For an ideal gas, the Mach number, Mach angle, shock- 
wave angle, and pressure ratio can be calculated at the leadin g 
edge, using the standard Rankine-Hugoniot shock relations 
and ut ilizing the free-stream Mach number and leading- 
edge deflection angle. With these flow parameters known, 
the only terms to be determined in equations (38), (39), 

(40), (43), and (44) are and^^ - These derivatives 

are easily obtained from the shock relations. (See ref. 9.) 


fdP\ (dP/da), 

\dS )„ ( d8/dcr)c 

where 

/ dP\ 2y Mg sin 2 a 

\d<r) t r 7+1 

and 

=1 ~ 2 ^?+x^ 0081 (<r--5)+sin 2(cr — 5) cot 2<r 


(Dl) 

(D2) 

(D3) 


By again using the standard forms of the . Ran k ine- 
Hugoniot shock relations, it is possible to transform equations 
(38), (39), and (43) (given in the analysis) into 


1 dP y tan f 

r j+i, 

^cos 2 /3 1 \ 

1 + 2 ! 

Vcos 2 f M<? sin 2 < 7 / 

K w dW sin 2 /3 

tan 2 f . 

s' 

Ceos 2 /3 1 \ 

tan 2 /S 

Vcos 2 f Mq 2 - sin 2 a) 


D 


2y Mr? sin 2 <r— (7— 1) 

7+1 


’] 


(D4) 


derivatives are most easily deterrriined by employing the 
temperature as a parameter. 18 Thus, 


and 


fdP\ (dP/dT ). 
\d5j' ( dS/dT X 

fdS\ _ (dS/dT), 
\daja (d(r/dT) a 


(D7) 

(D8) 


Differentiation of equation (18) yields 

pi T 0 (P-l-y 0 M 0 *)-M'TS (M. Jj+2y„ ~f) 


where from equation (15) 


T'+TJP 1 


~M*T, dy , 7 MS 


(D 9) 


dT ~\y,TMJ L : 

2 

dT' 2 1 

c e IT. ( ^ 

2 

y< 1 


) i 

7<-iJ 

and from equation (14) 

dy {yt 7) (7 1) | 

c 2 

6 e°i T ‘+ 1\ 

dT T„(y t 1) 1 

l 2 

T, e'lT—i) 


(DIO) 


(Dll) 


( d8 \ 

jjij can be evaluated from equation (20) 
7o-Mo 2 tan<r 1 




CP— l)* 


(S). 


for the surface pressure gradient, 


jB hr 7+1 

K„ 4 cos f 


1 tan 2 f 

tan 2 fi 

tan 2 f 1 /eos 5 /3 1 \ 

tan 5 (3 '2 Vcos 2 f Mo sin 2 a). 


(D5) 


for the shock-wave curvature, and 


1 + 




±0 

2 V 1 


cos 2 j3 . 1 

cos 5 f Mq 1 sin 2 <r 


) 


tan 2 f 
|_tan 2 (3 


1 / cos 2 /3 . 1 \ 

2 \cos 2 f ' Mo sin 2 cr/_ 


tan f 


tan j3 


(D6) 


for the surface-pressure-gradient ratio. These equations are 
similar to those given by Schaefer in reference 8 and require 
less work to compute than equations (38), (39), and (43). 

For calorically-imperfect-gas flows, the standard shock 
relations are obviously not applicable. In this case, the flow 
parameters at the leading edge can be evaluated using the 
oblique-shock-wave equations previously, presented. (See 
eqs. (14), (15), (18), (19), and (20).) Since the .primary 
variable in these equations is temperature, the required 
362 


yo-Mp 2 

P-1 



(D12) 


where from - equation (19) 



and from equation (8) 

(jr) a = T' \_ poT ° (jr)r p, l (D14) 

The procedure for calculating -i- 4yL> + and k is straight- 

jtXtP ^ V* JiX.lT 

forward with the aid of the preceding equations and may be 
summarized as follows: 

1. For any Mo, T 0 , 8^, and K a choose a value of T,/T 0 . 

i> The forms ol some of the derivatives presented have been somewhat simplified from the 
forms originally presented In reference 17. 
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(Fig. 0 or the ideal-gas relations will provide an initial 
estimate.') 

2. Calculate y c (eq. (14)) 

3. Calculated, (eq. (15)) 

4. Calculate P„ (eq. (18)) 

6. Calculate p 0 /p<, (eq. (8)) 

6. Calculate <r (eq. (19)) 

7. Calculate S„ (eq. (20)) 

If this value of 5» is not close enough to the desired value of 
<5.v, iterate, choosing a new TJT 0 . 

8. Calculate ^ (eq. (Dll)) 

9. Calculate (eq. (DlO)) 

10. Calculate (jiff) (eq. (D9)) 

11. Calculate (eq. (D14)) 

12. Calculate (eq- (D13)) 

13. Calculate (eq- (D12)) 

14. Calculate (jjjj (°q- (D7)) 

15. Calculate (eq. (40)) 

16. Calculate -gr- -j— (eq. (38)) 

17. Calculate i {/ (eq. (43)) 

18. Calculate (jfe) (eq. (D8)) 

19. Calculate (eq. (39)) 

-Ql to 

20. Calculate k (eq. (44)) 
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TABLE I.— FUNCTIONS FOR BLENDER-AIRFOIL METHOD 



7“ 

7< 

7a 

7“T« 


7“ 

7i 

To 

T 

“7a 

j(AfoJy) 



/(Aftf*) 

o(.McPh) 

/(AfoJjv) 







0 

L 000 

0 


1 

0 




0 



0 

.05 

L072 

.009901 


Him 

.009897 


5.2 

47.56 

.3487 

1.848 

44.61 

.3275 

.10 

L 148 

.01961 


L148 

.01959 


5.4 

5L13 

.3506 

L 347 

4a 33 


.15 

1.230 

.02912 



■ .02909 


5.6 

5182 

.3523 

L 347 

6L68 

.3309 

.20 

L316 

.03845 


L310 

.03841 


6.8 

58.66 


L 340 

55.20 

.3323 

.25 

1.406 

.01700 



.04755 


6.0 

62.62 

.3553 

L 346 

5a 90' 


.30 

1. 602 

.05656 

L 399 


.05649 


6.2 

66.73 

.3566 

L346 

62.79 

.3350 

.35 

LG04 

.06537 

L399 


.06524 


0.4 

7a 96 

.3*78 

L346 

06.77 


.40 

L 710 

.07393 

1.399 


.07380 


0.0 

75.33 

.3589 


.45 

1.823 

.0823* 

L 399 


.08220 


0.8 

79.83 

.3600 




.50 

L941 

.09058 

L399 


.09040 


7.0 

84.47 

.3610 




.55 

2.065 

.09863 

1.399 

2.064 

.09841 



sa 24 

.3618 


- 


.60 

2.195 

.1071 

L 399 

2.194 

.1062 


Brl! 

9114 

.3626 




.55 

2.332 

.1142 

1.399 

2 330 

. 1139 



9a 19 





.70 

2.474 

.1217 

L 399 

2 473 

.1213 



1014 





.75 

2.624 

.1290 

L 398 

2622 

.1286 



. ioa7 

.3647 




.80 

2.780 

.1362 

L39S 

2777 

.1356 



123.5 

.3661 




.85 

2.943 

.1431 

L 393 

2939 

.1425 



138.2 

.3673 




.90 

3.112 

.1499 

1.398 

3.108 

.1492 



153.8 

.3684 




.95 

a2S9 

.1565 


3.284 

.1557 



17a 2 

.3693 




LOO 

3.473 

.1630 

L 397 

3.466 

.1620 


i MM 

187.4 





LI 

3.862 

.1753 

L 396 

3.852 

.1740 


1L0 

205.4 

.3707 




L2 

4. 280 

.1869 

L 396 

4.266 

.1853 


11.5 

2313 

.3714 




L3 

4.728 

.1978 

L 395 

4.708 

1 .1958 


12.0 

2411 

.3719 




L4 

5.206 

.2064 

1.393 

5.179 

.2056 


12.5 

2617 

.3723 




L5 

5.715 

.2178 - 

L 392 

5.679 

.2147 


13.0 

286.1 

.3727 




L 0 

6.256 

.2269 

L 391 

6.207 

.2232 


13.6 

- 308.3 

.3731 




L 7 

6.827 

.2354 

L 389 

6.754 

.2309 


14.0 

33L4 

.3735 




L8 

7.431 

.2433 

1.388 

7.349 

.2381 


14.5 

355.4 

.3738 




L9 

8.066 

.2507 

L 386 


.2447 


16.0 

38a 2 

.3740 




2.0 , 

8,734 

.2577 

L 384 

8.605 

.2508 


' 16.0 

432.3 

.3745 




21 

9.411 

.2646 

L 382 

2 274 

.2564 


17.0 

487.7 

.3749 




2.2 

ia 17 

.2702 

1.380 

9.973 

.2616 


18.0 

54a 5 

.3752 




2.3 

10.93 

.2759 

1.378 

ia?o 

.2664 


19.0 

608.6 

.3756 




2.4 

1L 73 

.2812 

1.376 - 

1L46 

.2709 


20.0 

6712 

.3757 




2.5 

12.56 

.2862 

1.374 

12 24 

.2749 


22.0 

<815.3 

.3701 




2.0 

13.42 

.2908 

1.372 

13.06 

.2788 


24.0 

9018 

.3764 




2.7 

14.32 

.2951 

L 370 

12 90 

.2823 


26.0 

1138 

.3766 




Z8 

15.25 

.2992 

1.309 

14.77 

.2856 


28.0 

1319 

.3708 




2.9 

16.21 

.3030 

1-367 

15.68 

.2887 


80.0 

1614 

.3770 




ao 

17.21 

.3066 

L 365 

16.61 

.2916 


35.0 

2060 

.3772 




12 

10.30 

.3130 

L 302 

12 56 

.2968 ' 


4a 0 

2690 

.3774 




3.4 

2L52 

.3187 

L 360 

20. 64 

.3016 


4S.0 

3404 

.3776 




a« 

2a 88 

.3237 

L 358 

22.83 

.3057 


sao 

4202 

.3776 




as 

26.37 

.3280 

1.366 

25.16 

.3094 


6a o 

6050 

.3777 




4.0 

29.00 

.3322 

1.354 

27.60 

.3128 


7ao 

8233 

.3778 




4.2 

3L76 

.33*8 

L 352 

3a 17 

.31*8 


sao 

10750 

.3778 




4.4 

34.65 

.3389 

L361 

32.88 

.3186 


oao 

13610 

.3779 




4.6 

37.68 

.3418 

L350 
L 349 

3*. 68 

.3211 


. 100.0 

16S00 

.3770 




4.8 

4184 

.3443 

3a 63 

.3234 




.3780 




8.0 

44.14 

.3466 

L349 

41.70 

.3255 

_ 








TABLE II.— SURFACE PRESSURE GRADIENT, -g- ^ 


Mo 

Str 

- L6 

ao 

ZA 

10 

5.0 

ao 

ao 

10 JO 

15.0 

20 J& 

00 

0° 

2.817 

a 233 

1455 

5.784 

7.144 

a 5i9 

1L 29 

1107 

21.05 

2a 04 

a> 

2.0° 


a 530 


a 809 

a 767 

ia87 

15.69 

2a 99 

37.78 

*9.81 

CD 

*-0° 

a 616 

1014 


a 577 

1L 64 

ia2o 

2a 94 

35.02 

7101 

13X2 

CO 

7.5° 

1298 

1458 

a 896 

ia 24 

1141 

19.45 

32.30 

4a io 

109.3 

195.9 

CO 

iao° 

5.878 

1950 

7.849 

12.05 

17.46 

2112 

4L43 

6123 

145. 7 

20L6 

00 

15.0° 


a 185 

9.908 

15.96 

2100 

3106 

ea32 

9176 

21X1 

387.1 

00 

2a o° 


a765 

12.12 

ia©9 

3a 60 

4a 87 

7a 34 

12X2 

28X2 

50X7 

CO 

25.0° 



1156 

2a 94 

3a 80 

62.88 

9144 

14a 3 

33X2 

699.6 

CD 

3a o° 



ia25 

27.94 

42.54 

6a 91 

1083 

169. 7 

38X4 

682.8 

CO 

35.0° 




3a 89 

4& 83 

ea so 

12a 9 

m 5 

42X8 

75X6 

CO 

4a o° 





7L70 

8190 

139.3 

212.8 

471. 1 

83X8 

CO 

4a o° 









1727. 

1711. 

00 


TABLE m.— SURFACE-PRESSURE-GRADIENT RATIO, ^ 


X. Mo 

1.5 

ao 

XO 

4.0 - 

ao 

XO 

8 J) 

10.0 

15.0 

20.0 

00 

0° 


WWWB 


L000 

L000 

L 000 

L000 


L000 

1.000 

BSjJgiB 

zo° 


BlmB 

BV« 

L 000 

.9999 

.9998 

.9996 

.9993 

.9976 

.9918 


*.0° 



.9998 

.9994 

.9987 

.9978 

.9950 


.9773 

.9619 


7.5° 

L001 


.9995 

.9981 

.9961 

.9935 

.9804 

.9777 

.9560 

.9362 

.8825 

iao° 

.9978 

L001 

.9988 

.9959 

,9920 

.9870 

.9754 


.9309 

.9198 


15.0° 


LOOS 


.9895 

.9807 

.9715 

.9538 


.9168 

.9041 

.8846 

2 a o° 


L 025 

.9958 

.9824 

.9690 

.8578 

.9386 

.9253 

.9071 

.8991 

.8871 

2S.0° 



.9997 

.9783 

.9623 

.9493 

.9309 

.9194 

.9052 

.8994 

.8911 

3a o° 




.9825 

.9622 

.9182 

.9305 


.9087 

.9042 

.8979 

35.0° 



r 

a oi6 

.9787 

.9003 

.9407 


.9199 

.9158 


4a o° 





L 114 

1.028 

.9835 

.9672 

.9523 

.9472 


45,0° 









L737 

L325 

L 188 
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TABLE IY.— SHOCK- WAVE CURVATURE, K./K„ 


&N 

1.6 

2.0 

3.0 

4.0 

00 

6.0 

8.0 

10.0 

15J3 

200) 

as 

0° 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

a sooo 

2.0° 

.06597 

.04547 

. 04808 

.05767 

.06881 

.08080 

.1054 

.1305 

.1831 

.2537 

.8003 

8,0° 

.1951 

.1175 

.1212 

.1450 

.1728 

.2017 

.2508 

.3155 

.4357 

. 5264 

.8018 

7.8° 

.3698 

.1827 

.1827 

.2170 

.2566 

.2969 

.3733 

.4406 

.8683 

.6417 

.8042 

10.0° 

. 7988 

.2883 

.2447 

.2m 

.3354 

.3829 

.4676 

.5356 

.6460 

.7042 

.8077 

16.0° 


.4487 

.3704 

.4183 

.4737 

.5244 

.6044 

.6597 

’ .7345 

.7677 

.8181 

20, 0° 


9. 191 

.5051 

.6399 

.5894 

.6332 

.6970 

.7371 

.7863 

.8064 

.8352 

28.0° 



.6756 

.6634 

.6956 

.7271 

.7723 

.7997 

.8318 

.8445 

.8621 

3a 0° 



LOU 

.8222 

.8168 

.8294 

.8536 

.8696 

.8885 

.8961 

.9064 

38.0° 




LX68 

L 021 

.9890 

.9774 

.9779 

.9818 

.9840 

.9876 

4a 0° 





2.040 

L402 

L 203 

L 240 

L 202 

L 191 

L 179 

45,0° 









7.668 

3.829 

2.811 


TABLE V.— SHOCK- WAVE-CURVATURE RATIO, k 


Mo 
&N N. 

1.5 

10 

3.0 

4.0 

5.0 

6.0 

8.0 

10.0 

> 

15.0 

20.0 

CD 

0° 

L 000 

1.000 

L000 

L 000 

LOGO - 

L 000 

L 000 

L000 

LOCO 

L000 

1.072 

10° 

.9997 

.9998 

L 000 

L 000 

1.001 

L001* 

1.002 

L 003 

1.007 

1.012 

L072 

8.0° 

.9984 

.9992 

L 001 

L 002 

L 004 

L 006 

L011 

L 017 

L 030 

L 040 

L072 

7.5° 

.9970 

.9983 

1.002 

L 005 

L 009 

L 013 

L 021 

L 029 

LOU 

L 053 

1. 071 

10.0° 

1.002 

■ 996S 

1.003 

L 009 

L 014 

L 020 

L 030 

L038 

L 052 

L059 

1.071 

16.0° 


.9912 

1.005 

L 015 

Uffl 

L 031 

1.042 

L049 

L059 

L 063 

1.069 

2 a 0° 


.9795 

L 005 

L 019 

L029 

L 037 

L047 

L 053 

L 060 

L063 

L067 

28.0° 



1.000 

L019 

1.031 

L039 

1.048 

L053 

L058 

L 060 

1.063 

3a 0° 



.9830 

L013 

1.027 

L 035 

L044 

L 048 

1.053 

L055 

L057 

38.0° 




.9909 

L 012 

L 023 

L034 

L 039 

L044 

1.048 

L048 

4a 0° 





.9887 

.9879 

1.008 

L 016 

L023 

1.025 

L028 

48.0° 









.8861 

.924 2 

.9476 


TABLE VI.— SURFACE PRESSURE GRADIENT, SURFACE- 
PRESSURE-GRADIENT RATIO, SHOCK-WAVE CURVA- 
TURE, AND SHOCK-WAVE-CURVATURE RATIO FOR A 
CALORICALLY IMPERFECT, DIATOMIC GAS 

[To-«0° Hankino] 


M 0 

6° 

T,IT 0 

1 dP 
K w dW 

* 

K,IK U 

K 

3 

30.48 

2.01200 

18.36 

L 026 

L021 

0.9828 

6 

8.042 

1.32750 

16.04 

.9954 

.2732 

L010 

6 

15. 30 

L 75357 

24.38 

.9796 

.4773 

L025 

6 

2a 48 

2.13800 

30.94 . 

.9075 

.5866 

L031 

5 

2 a oi 

2.62500 

37.39 

.9577 

.6908 

L034 

6 

3L09 

3.12700 

42.63 

.9549 

.7986 

L031 

6 

37.30 

3. 81250 

49.60 

.9725 

L041 

L015 

6 

41.01 

4.30000 

6a79 

L 042 

L 571 

.9827 

10 

10.01 

2.10800 

63.86 

.0621 

.5265 

1 . (M0 

10 

10.38 

2.17000 

sa 07 

.9598 

.5373 

L041 

10 

20.62 

4.40000 

123.6 

.0112 

.7132 

L 082 

10 

27.61 

a 45000 

166.2 

.8932 

.7922 

L 065 

10 

3L27 

7.64000 

169.8 

.8886 

.8363 

L 064 

10 

4a is 

•la 87600 

195.0 

.8969 

L018 

L 060 

10 

42.98 

•12.00000 

204.1 

.9128 

1.151 

L039 

15 

5.230 

L 77500 

7a 83 

.9750 

.4464 

L032 

16 

ia49 

3.25000 

me 

.9283 

.6363 

L 067 

15 

16.31 

6.60000 

229.5 

.8936 

.7250 

L 070 

15 

22.28 

8.70000 

30L4 

.8733 

1 7834 

L074 

15 

28.65 . 

* 12.75000 

364.5 

.8630 

.8391 

L074 

20 

6.122 

2. 15000 

132.4 

.9594 

.5242 

L 042 

20 

iae7 

4.70000 

274.2 

.9027 

.6923 

L068 

20 

ias3 

8.70000 

417.0 

.8706 

.7648 

L 077 

20 

2a 98 

•12.60000 

514.4 

.8579 

.8011 

L079 

20 

32.61 

•26.80000 

705.0 

.8490 

.8914 

L 076 

20 

39. 15 

•36.00000 

77L1 

.8571 

.9861 

L066 


» Values ot the temperature greater than 5,000° E. downstream of the shook wave. 








